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Field of the Invention 

[OOOl] The present invention relates generally to improvements in radio 
positioning systems and their methods of operation and, more particularly, to methods 
and apparatus for simplifying the acquisition of data required in such systems. 

Background of the Invention 

[0002] EP-A-0 303 371, the contents of which are hereby incorporated by 
reference, describes a radio navigation and tracking system which makes use of 
independent radio transmitters set up for other purposes. The signals from each 
transmitter, taken individually, are received by two receiving stations, one at a fixed 
and known location, and the other mounted on the mobile object whose position is to 
be determined. A representation of the signals received at one receiving station is sent 
via a link to a processor at the other receiving station, where the received signals are 
compared to find their phase differences or time delays. Three such measurements, 
made on three widely spaced independent transmitters, are sufficient to determine the 
position of the mobile receiver in two dimensions, i.e. its position on the ground. The 
phase or time offset between the master oscillators in the two receivers is also 
determined. 

[0003] "CURSOR", as the system described in EP-A-0 303 371 is known, is a 
radio positioning system which can use the signals radiated by existing non- 
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synchronised radio transmitters to locate the position of a portable receiver. Unlike 
some other systems which use the temporal coherence properties of networks of 
purpose-built synchronised transmitters, CURSOR makes use of the spatial coherence 
of the signals transmitted by single transmitters. In a further development (see EP-A-0 
880 712 & WO-A-99/21028), the technology has been applied to find the position of a 
mobile phone handset in a GSM or other digital telephone system, and these are 
examples of an 'Enhanced Observed Time Difference' (E-OTD) method using the 
down-link signals radiated by the network of Base Transceiver Stations (BTS) of the 
telephone system. 

[0004] In the digital mobile telephone application described in EP-A-0 880 712, 
the contents of which are hereby incorporated by reference, the signals from each BTS 
within range of the handset are received both by the handset itself and by a fixed 
nearby receiver, the Location Measurement Unit (LMU), whose position is accurately 
known. Representations of the received signals are passed to a Mobile Location 
Centre (MLC) where they are compared in order to find the time difference between 
them. Figure 1 shows the geometry of a standard two-dimensional system. The origin 
of Cartesian co-ordinates x and y is centred on the LMU positioned at O. The 
orientation of the axes is immaterial, but may conveniently be set so that the y axis lies 
along the north-south local map grid. The handset, R, is at vector position r with 
respect to the LMU position O. A BTS, A, is shown at vector position a. 
[0005] Consider first the signals from BTS A. The time difference, Ah, measured 
between the signals received at R and O is given by 

At a =(\r-a\-\a\)/o + e, 

where v is the speed of the radio waves, s is the clock time offset between the clocks 
in the receivers at R and O, and the vertical bars each side of vector quantities denote 
that it is the magnitude of the vectors which are used in the equation. The value of e 
represents the synchronisation error between the measurements made by the two 
receivers. Similarly, may be written for two other BTSs (B and C) at vector positions 
b and c (not shown): 

Af b =(|r-b|-|b|)/u + f, 
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A/ c =(|r-c|- |c|)/l> + €. 



(1) 



[0006] The values of At a , At b , At c , are measured by the methods disclosed in EP- 
A-0 880 712 and the values of a, b, c, and v are known. Hence the equations (1) can 
be solved to find the position of the handset, r, together with the value of s . 
[0007] In WO-A-99/21028, the contents of which are hereby incorporated by 
reference, it is described how these same time offsets can be measured using locally- 
created templates in a GSM telephone system as follows. Suppose that the handset R 
has recorded a short burst of the GSM signals from BTS A. Contained within that 
recording is the framing structure, synchronisation bursts and other 'given' data (or 
predetermined values) which are a constant feature of those transmissions. The 
processor within the handset can create a matching template, based on the known 
structure of the network signals. Received signals can then be matched by the locally- 
generated template. When the template finds a match, the correlation peak at the 
position of best match corresponds to the time offset between the received signals and 
the local clock inside the handset. For the signals radiated by BTS A this measured 
time offset, Af a i, is given by 

Af a i =(|r - a|)/u + S\ , 

where Oq is the time offset of the BTS transmissions and €\ is the time offset of the 
handset's internal clock, both relative to an imaginary universal 'absolute' clock. The 
signals from BTSs B and C may also be measured in the same way, giving 

A*bi =(|r-b|)/u+ <2b+ s u 

and 

Afci =(|r - c\)/o + s x . (2) 
[0008] The same measurements can also be made by the LMU, giving 
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A*a2 =(|a|)/u + «a+ ft , 



A/ b 2 =(|b|)/y + £2 , 

and 

Afe =(|c|)/u + Oc+ ft , (3) 

where £2 is the time offset of the LMU's internal clock relative to the same imaginary 
universal absolute clock. Subtracting equations 3 from equations 2 gives 

At a = Afai - A^2= (|r - a| - |a|)/u + s , 

A* b = Aki - A^ b 2= (|r - b| - |b|)/t> + £ , 

and 

At c = A^ci - Af C 2= (|r - c| - |c|)/u + e , (4) 

where £ = £1 - ft. It will be noted that equations 4 are just like equations 1, and can be 
solved in the same way to find the position of the handset, r, and the value of s . 
[0009] It will be apparent that the CURSOR method as described above, in 
common with all other methods which use the signals from non-synchronised 
transmitters, requires a network of LMUs to be set up within the coverage area of the 
telephone system. These units act as reference points at which the unsynchronised 
signals radiated by the BTSs are measured for comparison with the same signals 
received by a handset. Each position measurement requires a match to be made 
between the signals received by the handset from a number of nearby BTSs, and 
signals received by an LMU from the same set of BTSs. In practice, it is often 
difficult to find a match using just one LMU, especially if the LMU network is sparse, 
since the handset may receive signals from BTSs not received by the LMU, and vice- 
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versa. It is therefore necessary to combine the measurements from two or more 
LMUs. However, each new LMU brought into the calculation adds a further unknown 
clock time offset (s 2 , £3 etc.), each of which therefore requires an additional BTS 
measurement to provide the extra equation needed to solve for all the unknown 
quantities. 

[OOIO] One solution to this problem is presented in WO-A-99/21028 where it is 
shown how the LMU network can be synchronised. Referring to Figure 2, suppose 
that an adjacent pair of LMUs, Ui and U2, can see a common BTS. The positions of 
the LMUs and the BTS are all known, so a single measurement of the BTS signals by 
each LMU is sufficient to determine the clock time offset between the LMUs. For 
example, suppose that the distance from Ui to the BTS is s\, and the distance from U2 
to the BTS is j 2 - Ui measures time offset A*i and U 2 measures At 2 , given by 

A*i = s\lu + a + £21 > 

A*2 = s 2 lo + a + £22 , (5) 

where a is the time offset of the BTS transmissions, and s 2 \ and £22 are the time 
offsets of the LMU internal clocks in Ui and U2 respectively. Subtracting the second 
equation from the first yields 

&i\ - £22 ~ A*i - At 2 + s\lv - s 2 lv , (6) 

which is the relative time offset of the clock in Uj with respect to that in U2. This 
process may be repeated for a second pair of LMUs, say U2 and U3, and another BTS 
whose signals can be received by both members of this second pair of LMUs. In this 
way a synchronisation map may be calculated, which provides the clock offsets of all 
the LMU internal clocks relative to one of them adopted as a master C LMU network 
clock time'. Having established the LMU synchronisation map in this fashion, a 
CURSOR position measurement can then include any number of LMUs without the 
penalty of adding an extra unknown time offset for every LMU, since the relative 
LMU time offsets are known. 
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[OOl 1] The receivers discussed in the preceding paragraphs make measurements of 
time offsets. More generally, receivers can measure time offsets, phase offsets (which 
can be converted into time offsets with a modulo 360° ambiguity), frequency offsets or 
rates of change of frequency offsets. Though these measurements are of different 
quantities, the present invention is applied usefully to each of them as, when combined 
with similar measurements made by a second receiver, they can independently provide 
positional information. Positioning systems making use of these measurements are 
discussed in a related patent application (our reference MJB06427WO) filed 
simultaneously herewith. 

Summary of the Invention 

[0012] The present invention teaches how the same advantages of an (effectively) 
synchronised LMU network may be obtained by setting up one or more 'virtual 
LMUs' in the network which act as interface nodes for the real LMUs. 
[0013] According to a first aspect of the invention, there is provided a method of 
generating a list of offsets in time, phase, frequency, or derivatives thereof, or their 
equivalents expressed as offsets in distance or derivatives thereof, of a plurality of 
transmission source signals, corresponding to a given location, relative to a common 
reference, the method comprising 

(a) acquiring data from two or more receivers, the positions of which may 
be known or determined, the data from a receiver comprising offsets in 
time, phase, frequency, or derivatives thereof respectively of signals 
received from the transmission sources relative to a reference source in 
each receiver or to each other; and 

(b) combining the acquired data and calculating the list of offsets 
corresponding to the given location relative to the common reference. 

[0014] In practice the offsets from the list can be used in place of offsets obtained 
directly from the receiver or receivers. 

[0015] The relative offsets in time, phase, frequency, or derivatives thereof, with 
respect to each other or with respect to a reference source, of the signals received by a 
first receiver from a plurality of the transmission sources may be represented by 
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corresponding offsets or differences in the distances between the transmission sources 
and the first or second receivers. 

[0016] The invention also includes apparatus using the method above, the 
apparatus comprising 

(a) means for acquiring data from two or more receivers, the positions of 
which may be known or determined, the data from a receiver 
comprising offsets in time, phase, frequency, or derivatives thereof 
respectively of signals received from the transmission sources relative 
to a reference source in each receiver or to each other; and 

(b) means for combining the acquired data and calculating the list of 
offsets corresponding to the given location relative to the common 
reference. 

[0017] In a method using techniques similar to or as described in EP-A-0 880 712, 
in place of offsets in time, phase, frequency, or derivatives thereof, or their equivalents 
expressed as offsets in distance or derivatives thereof, data representative of the 
received signals may be used, from which the offsets of signals received from the 
transmission sources relative to the reference source may be determined. 
[0018] Therefore, the invention also includes a method of generating a list of 
offsets in time, phase, frequency, or derivatives thereof, or their equivalents expressed 
as offsets in distance or derivatives thereof, of a plurality of transmission source 
signals, corresponding to a given location, relative to a common reference, the method 
comprising 

(a) acquiring data from two or more receivers, the positions of which may 
be known or determined, the data from a receiver being representative 
of the received signals; 

(b) determining from the acquired data the offsets in time, phase, 
frequency, or derivatives thereof respectively of signals received from 
the transmission sources relative to a reference source or to each other; 
and 

(c) combining the offsets so determined and calculating the list of offsets 
corresponding to the given location relative to the common reference. 

[0019] The invention also includes apparatus for carrying out the method 
described immediately above, the apparatus comprising 
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(a) means for acquiring data from two or more receivers, the positions of 
which may be known or determined, the data from a receiver being 
representative of the received signals; 

(b) means for determining from the acquired data the offsets in time, phase, 
frequency, or derivatives thereof respectively of signals received from 
the transmission sources relative to a reference source or to each other; 
and 

(c) means for combining the offsets so determined and calculating the list 
of offsets corresponding to the given location relative to the common 
reference. 

Brief Description of the Drawings 

[0020] Figure 1 shows the geometry of a CURSOR system as described in EP-A-0 
880 712; 

[002 1] Figure 2 shows adjacent LMUs making measurements of a common BTS; 
[0022] Figure 3 shows a network of real and virtual LMUs in a system of the 
invention; 

[0023] Figure 4 illustrates a similar, simplified network; 

[0024] Figure 5 shows the positions of LMU sites and BTS sites in a real network; 
[0025] Figure 6 shows the same network with the addition of a Virtual LMU; 
[0026] Figure 7 shows the same network with the real LMUs replaced by the 
Virtual LMU; 

[0027] Figure 8 illustrates, by way of a flow chart, the processing which takes 
place inside an LMU; 

[0028] Figure 9 illustrates, by way of a flow chart, the processing which takes 
place inside the VLMU; 

[0029] Figure 10 shows a list of timing offsets generated in a real system such as 
that shown in the example of Figures 5 to 9; and 

[0030] Figure 1 1 shows a table of timing errors associated with LMUs and a 
virtual LMU in this example. 
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Detailed Description of the Exemplary Embodiments 

[0031] A radio positioning method and system including either of the methods 
and apparatus defined above also form part of the present invention. 
[0032] The invention also includes apparatus (a 'virtual LMU') for carrying out 
either or both of these methods. The apparatus may include a computer (located 
anywhere convenient) and programmed to carry out the required process. Although 
the following description of a particular application of the invention concerns signals 
in a digital telephone network, it will be apparent that the invention is by no means 
restricted to this application but may be applied to any network of one or more 
transmitters, synchronised or unsynchronised, set up for any purpose. 
[0033] A virtual LMU includes a computer process which may run on any 
computer platform able to obtain data from real LMUs. Accordingly, a further aspect 
of the invention includes a method of calculating and maintaining a list of offsets in 
time, phase, frequency, or derivatives thereof, or their equivalents expressed as offsets 
in distance or derivatives thereof, of a plurality of transmission source signals, 
corresponding to a given location, relative to a common reference. 
[0034] It is assumed that the network of BTSs is unsynchronised in that the 
transmission time offsets of the BTS signals bear no constant or known relationship to 
each other, but that nevertheless the BTS oscillators are quite stable 9 so that their 
instantaneous frequencies change only slowly with time. In these circumstances, it is 
possible to predict the currently received offset in time, phase, frequency, or 
derivatives thereof of the signals from a given BTS by a given real LMU from 
sufficiently-recent historical data. The real LMUs in the network make measurements 
of all the BTSs they can detect in a cyclic fashion, repeating the cycle every few 
seconds. They maintain these measurements in a stack, replacing the oldest 
measurements with the most recent. A linear or low-order polynomial fit to the 
measurements therefore provides a predictor for extrapolation into the near future, or 
for interpolation in the recent past. Assume that the BTS oscillators are sufficiently 
stable that reliable predictions can be made over a period of, say, ten minutes. Then, 
every few minutes, the virtual LMU (VLMU) contacts each real LMU and receives its 
predictors for the received offsets of the signals from all the BTSs in its measurement 
set. It is likely that many of the BTSs will have been measured by more than one 
LMU, so the VLMU analyses the complete data set from all the real LMUs using well- 
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known methods to determine both (a) the best values of the real LMU internal clock 
offsets in time, phase, frequency, or derivatives thereof, and hence (b) the received 
offsets in time, phase, frequency, or derivatives thereof of the signals from all the 
BTSs which would have been measured by a real LMU located at the assumed 
position of the VLMU and capable of receiving signals from every BTS. 
[0035] In the above description of the function of the VLMU, it should be 
understood that any or all of the LMUs could be replaced by other receivers, not 
necessarily fixed or at known positions, which have not been set up specifically as 
LMUs. For example, the data from a number of handsets could be used to determine 
frequency offsets if the handsets were stationary. Furthermore, it is shown in a related 
patent application (our reference MJB06427WO) filed simultaneously herewith, how 
the positions and velocities of handsets can be determined without the need for any 
LMUs at all. 

[0036] Particular advantages of using a VLMU in a network include the following: 

(a) a full match can be guaranteed between a handset's measurements and a single 
(virtual) LMU; 

(b) the VLMU process minimises the timing errors in individual LMU measurements; 

(c) the VLMU list is immediately available to the position calculation process, 
increasing the speed of computation; 

(d) when combined with the ideas described in a related patent application (our 
reference MJB06427WO) filed simultaneously herewith, a list of receive time 
offsets can be created in a network of BTSs where there are few if any real LMUs; 

(e) the VLMU provides, in effect, a synchronisation map of the real LMU network 
which network may then be used to monitor the BTS network and, in particular, to 
determine locations of newly-installed BTSs. 

[0037] One example of a method and apparatus according to the present invention 
will now be described with reference to the accompanying drawings, in which: 
Figure 1 shows the geometry of a CURSOR system as described in EP-A-0 880 712; 
Figure 2 shows adjacent LMUs making measurements of a common BTS; 
Figure 3 shows a network of real and virtual LMUs in a system of the invention; 
Figure 4 illustrates a similar, simplified network; 

Figure 5 shows the positions of LMU sites and BTS sites in a real network; 
Figure 6 shows the same network with the addition of a Virtual LMU; 
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Figure 7 shows the same network with the real LMUs replaced by the Virtual LMU; 
Figure 8 illustrates, by way of a flow chart, the processing which takes place inside an 
LMU; 

Figure 9 illustrates, by way of a flow chart, the processing which takes place inside the 
VLMU; 

Figure 10 shows a list of timing offsets generated in a real system such as that shown 
in the example of Figures 5 to 9; and 

Figure 1 1 shows a table of timing errors associated with LMUs and a virtual LMU in 
this example. 

[0038] By way of example, and with reference to Figure 3, the function of a virtual 
LMU is now described which determines timing offsets. 

[0039] Consider a network of N real LMUs and M BTSs which includes a virtual 
LMU (VLMU). The position of the n th LMU, U n? is represented by the vector u n and 
the position of the m th BTS, B m , is represented by the vector b m , both vectors being 
with respect to the same origin. Signals radiated by BTS m will be received by LMU 
n after a time lag, and the measurement of this time lag, At nm , is given by 

A* nm = |u n -b m |/0 + ft + «m ± ^nm , ( 7 ) 

where ft is the clock time offset of LMU n, a m is the transmission time offset of BTS 
m, both with respect to an imaginary universal 'absolute* clock, and a nm is an estimate 
of the error in the measurement of Ar nm . The assumed position of the virtual LMU, V, 
is represented by the vector v. If the VLMU were able to receive the signals directly 
from BTS m without error then it would measure a receive time offset, fi m9 relative to 
the imaginary universal absolute clock, given by 

An = iv-b m |/u + a m . (8) 

[0040] Substituting for Om in equation (7) using the value deduced from equation 
(8) gives 

Afnm = |u n -b m |/u + ft + fim - |V-b m |/u ± (W < 9 ) 
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[0041] Over the entire network of N LMUs, all of the M BTSs are visible. Each 
individual LMU, however, will only see a few of them, but as long as there is 
significant overlap of visibility, it is possible to take the set of all At values and solve 
for values of and fi m . Hence the VLMU can calculate timings for any BTS as if the 
network of LMUs were synchronised, or as if just one LMU (the VLMU, V) was all 
that was needed to cover the entire BTS network. 

[0042] To illustrate this further, a simplified problem is shown and solved below 
using N=2 and M = 4, i.e. a network consisting of just 2 LMUs monitoring 4 BTSs 
(see Figure 4). For simplicity it is chosen that si = 0. This is allowed as the 'absolute' 
clock time is completely arbitrary and may, for example, be measured by the internal 
clock of LMU number one. (Note, however, that this choice does introduce an 
asymmetry into the solution in that the error associated with BTSs three and four are 
not the same.) The first LMU (U,) can receive signals from BTS numbers one, two 
and three, but cannot receive the signals from BTS number four. The second LMU 
(U 2 ) can receive the signals from BTSs one, two and four, but cannot see BTS number 
three. The equations can be written in matrix form as 



1 0 0 0 0 

0 10 0 0 

0 0 10 0 

1 0 0 0 1 
0 10 0 1 
0 0 0 1 1 



At l2 - 



lu.-b.l-lv-b.l 

U 

|»i -b 2 |-|v-b 2 



v 



At v 



At, 



k -b, - v-b, 



lu.-b.j-lv-M 



At 22 -- 



lu, -b,|- v-b 2 



v 



At, 



■b 4 - v-b 4 



+ Z 



(10) 



or equivalently as A.x = b + Z, where Z is an unknown vector of the actual errors on 
each measurement. 

[0043] The standard technique known as 'least squares' postulates that the 
estimate for x which minimises Z is given by 

x = (A T WA)" 1 A T Wb, < n ) 
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where the symbol A T indicates the transpose of the matrix A, and the matrix W is 
defined by 



( 1 



w = 



' 11 
0 

0 

0 

0 

0 



0 

1 



0 
0 
0 
0 



0 
0 
0 



0 
0 
0 

1 

2 

0 
0 



0 
0 
0 
0 

1 



0 
0 
0 
0 
0 

1 



(12) 



[0044] This particular example can be solved explicitly. For simplicity, it is 
assumed that all the values of a nm are the same, and equal to a This gives the result 



fix = ( 3Ai + D X2 + D 2 \ -D 22 )/4 
^ 2 = (3A 2 + Ai + A2-A1 )/4 

fi 4 = ( 2D 24 + A1+A2- D 2 i - D 22 )/2 

S2=(D 2 l+£>22-Al-A2)/2 



± 0.87a ? 
± 0.87a 5 
± 1.00a, 
± 1.41a, 
± 1.00a, 



(13) 



where 



lu -b J- v-b, 



v 



(14) 



[0045] Note that even in this simple case, when a BTS is seen by both LMUs, the 
errors in the calculated offsets are less than those in each of the measurements 
themselves. This is an important advantage of the virtual LMU method. 
[0046] The LMUs may also contain other means of synchronisation. For example, 
each real LMU could be connected to a GPS or other timing reference receiver, G, 
which serves to provide the common timing reference. In this case, the LMU network 
may be considered as synchronised already to this common timing reference (say GPS 
standard time), and then the VLMU need not solve for the individual values of s as 
these are already known. One advantage of using other means of synchronisation is 
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that there is no longer a requirement for overlap of visibility of BTSs between adjacent 
LMUs. If each BTS site also carried an LMU, then that LMU would need only to be 
able to receive the (very strong) signals from its co-sited BTS transmitter(s), thereby 
simplifying the installation of the LMU antenna. 

[0047] The VLMU mode of operation described above may be called the 'pull 
mode' as it requires the VLMU to instigate data transfer to itself from every real 
LMU. It is also possible to have each real LMU continuously check the difference 
between its own prediction of the receive time offset from every BTS using the values 
calculated from the set of predictors last sent to the VLMU and the actual measured 
values. When any one of these differences exceeds a given value, the LMU can send 
its new predictor set to the VLMU. This mode of operation may be called the 'push 
mode'. The particular mode appropriate for a real system depends, amongst other 
things, on the stability of the BTS network. 

[0048] An example of a prototype system constructed in accordance with the 
invention will now be described, in which timing measurements made by a number of 
LMUs are combined to create a list of receive time offsets for all the GSM BTSs in 
and around Cambridge, UK, as if they were observed by a single "Virtual" LMU. 
[0049] The positions of nine LMUs, U, - U 9 (shown as filled circles) and twenty 
three Cambridge area BTSs, B, - B 23 (shown by unfilled squares) are plotted in Figure 
5 on an Ordinance Survey (OS) grid. Each LMU comprises (a) hardware including an 
internal clock, a GSM radio, a computer, and a telephone connection, and (b) software 
including a program to compile a list of receive time offsets. 

[0050] Figure 8 illustrates, by way of a flow chart, the main elements of the LMU 
software residing in each LMU, U, - U 9 . Every few seconds the "scan cycle" program 
is entered at step Al . The program proceeds, in steps A2, A3, A6 and A7, to tune the 
GSM radio to each GSM transmission channel in turn and scan for BCCH signals. In 
the event that a BCCH signal is detected in step A3, the program calculates the receive 
time offset of the signal relative to its internal clock (in step A4) and updates its list of 
time offsets accordingly in step A5. The BCCH is also decoded to produce the ID of 
the BTS B n from which the signal is received. Figure 5 also shows lines between 
LMUs and the BTSs detected by them as a result of the scan cycle routine. 
[0051] Figure 6 shows the positions of the nine LMUs, the twenty three BTSs 
monitored by those LMUs and a Virtual LMU V (plotted as an unfilled circle). The 
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Virtual LMU V comprises (a) hardware including a telephone connection to each real 
LMU and a computer, and (b) software including a program to compile a list of 
"virtual" receive time offsets. 

[0052] Figure 9 illustrates, by way of a flow chart, the main elements of the 
VLMU software. Every 4 minutes the "VLMU update" program is entered at step VI . 
The program proceeds, in steps V2, V3, V5 and V6, to connect to each of the real 
LMUs via the telephone connections. If the connection to a given LMU Ui - U9 
succeeds, the program, in step V4, extracts that LMU's list of receive time offsets. On 
completion of the connection cycle, the program (in step V7) combines the data to 
produce a list of receive time offsets for all the BTSs monitored by the LMU network. 
An example of this list is shown in Figure 10 where part of a table of receive time 
offsets generated by the VLMU is reproduced. The first column shows the BTS 
identifier and the figures in the table are in units of 1.85 microseconds. Figure 6 also 
shows lines between the VLMU V and each real LMU Ui - U 9 from which the VLMU 
software is able to extract a list of time offsets. 

[0053] Once the "VLMU update" has been completed, the combined list of 
timings generated is equivalent to that which would have been observed had there 
been a single (real) LMU at the VLMU's location making timing measurements of 
every BTS in the network. Figure 7 illustrates this equivalence by showing the virtual 
monitoring of each BTS's timing measurements by the Virtual LMU (shown as lines 
between the VLMU V and the BTSs Bi - B 23 )- 

[0054] Each timing measurement has an associated error which, in most cases, is 
smaller than the errors in the timing measurements made by the individual real LMUs. 
This is illustrated in the table of Figure 11, which shows part of the list produced by 
the VLMU during operation. The first column shows the identifier of each BTS. The 
columns headed 1 to 9 each relate to a particular one of the nine real LMUs. The 
figures in the table are errors in the timing measurements made by the real LMUs of 
the signals from the corresponding BTSs. An empty cell indicates that the LMU is 
unable to receive a BTS's signal. The column headed VLMU shows the result of 
combining the measurements using the method outlined above. 

[0055] The timings in the table of Figure 11 are in units of 1.85 microseconds. 
Note that the VLMU timing errors are generally smaller than those estimated for real 
LMU timings, confirming an advantage of the VLMU method in reducing errors. 
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What Is Claimed Is 

1. A method of generating a list of offsets in time, phase, frequency, or 
derivatives thereof, or their equivalents expressed as offsets in distance or derivatives 
thereof, of a plurality of transmission source signals, received at a given location, 
relative to a common reference, the method comprising 

(a) acquiring data from plural receivers, the positions of which may 
be known or determined, the data from a receiver comprising offsets in time, 
phase, frequency, or derivatives thereof, respectively of signals received from 
the transmission sources relative to a reference source in each receiver or to 
each other; and 

(b) combining the acquired data and calculating the list of offsets relative 
to the common reference. 

2. A method of generating a list of offsets in time, phase, frequency, or 
derivatives thereof, or their equivalents expressed as offsets in distance or derivatives 
thereof, of a plurality of transmission source signals, received at a given location, 
relative to a common reference, the method comprising 

(a) acquiring data from plural receivers, the positions of which may 
be known or determined, the data from a receiver being representative of the 
received signals; 

(b) determining from the acquired data the offsets in time, phase, 
frequency, or derivatives thereof, respectively of signals received from the 
transmission sources relative to a reference source or to each other; and 

(c) combining the offsets so determined and calculating the list of offsets 
relative to the common reference. 

3. A radio positioning method for determining the position of one or more 
receivers the positions of which are unknown, which method includes the method of 
claim 1 or claim 2. 

4. A radio positioning method according to claim 3, wherein the common 
reference comprises an external reference. 
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5. A radio positioning method according to claim 4, wherein the common 
reference comprises a GPS signal. 

6. A radio positioning method according to claim 3, wherein the step of acquiring 
data from said plural receivers includes instigating acquisition of said data from a 
common location. 

7. A radio positioning method according to claim 3, wherein the step of acquiring 
data from said plural receivers includes instigating acquisition of said data from each 
said receiver at times determined by each said receiver. 

8. Apparatus for generating a list of offsets in time, phase, frequency, or 
derivatives thereof, or their equivalents expressed as offsets in distance or derivatives 
thereof, of a plurality of transmission source signals, received at a given location, 
relative to a common reference, the apparatus comprising 

(a) means for acquiring data from plural receivers, the positions of 
which may be known or determined, the data from a receiver comprising 
offsets in time, phase, frequency, or derivatives thereof, respectively of signals 
received from the transmission sources relative to a reference source in each 
receiver or to each other; and 

(b) means for combining the acquired data and calculating the list of 
offsets relative to the common reference. 

9. Apparatus for generating a list of offsets in time, phase, frequency, or 
derivatives thereof, or their equivalents expressed as offsets in distance or derivatives 
thereof, of a plurality of transmission source signals, received at a given location, 
relative to a common reference, the apparatus comprising 

(a) means for acquiring data from plural receivers, the positions of 
which may be known or determined, the data from a receiver being 
representative of the received signals; 

(b) means for determining from the acquired data the offsets in time, phase, 
frequency, or derivatives thereof, respectively of signals received from the 
transmission sources relative to a reference source or to each other; and 



17 



(c) means for combining the offsets so determined and calculating the list 
of offsets relative to the common reference. 

10. A radio positioning system including apparatus according to claim 8 or to 
claim 9. 

11. A radio positioning system according to claim 10, wherein the common 
reference comprises a reference external to said receivers. 

12. A radio positioning system according to claim 11, wherein the common 
reference comprises a GPS signal. 

13. A radio positioning system according to claim 10, wherein the means for 
acquiring data from said plural receivers includes a computer system arranged to 
instigate the transfer of said data from said plural receivers to said computer system at 
times determined by said computer system. 

14. A radio positioning system according to claim 10, wherein the means for 
acquiring data from said plural receivers includes a computer system, and including 
means for instigating said acquisition of data from each said receiver at times 
determined by each said receiver. 

15. A digital telephone network, including a radio positioning system according to 
claim 10. 
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